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Abstract 

Background: Recent genome wide association studies (GWAS) and previous positional linkage studies have 
identified more than 50 single nucleotide polymorphisms (SNPs) associated with obesity, mostly in Europeans. 
We aimed to assess the contribution of some of these SNPs to obesity risk and to the variation of related metabolic 
traits, in Mexican children. 

Methods: The association of six European obesity-related SNPs in or near FTO, NPC1, ENPP1, NEGRI, GNPDA2 and 
MC4R genes with risk of obesity was tested in 1,463 school-aged Mexican children (W coses = 514; /V confro , s = 949). We 
also assessed effects of these SNPs on the variation of body mass index (BMI), fasting serum insulin levels, fasting 
plasma glucose levels, total cholesterol and triglyceride levels, in a subset of 1,171 nonobese Mexican children. 

Results: We found a significant effect of GNPDA2 rs1 0938397 on risk of obesity (odds ratio [OR] = 1.30; P= 1.34 x 10~ 3 ). 
Furthermore, we found nominal associations between obesity risk or BMI variation and the following SNPs: ENPP1 
rs7754561 , MC4R rsl 778231 3 and NEGRI rs281 5752. Importantly, the at-risk alleles of both MG4R rs1 778231 3 and NPG1 
rsl 805081 showed significant effect on increased fasting glucose levels ((3 = 0.36 mmol/L; P = 1.47 X 10~ 3 ) and 
decreased fasting serum insulin levels ((3 = -0.10 uU/mL; P= 1.21 x 10~ 3 ), respectively. 

Conclusion: Our present results suggest that some obesity-associated SNPs previously reported in Europeans also 
associate with risk of obesity, or metabolic quantitative traits, in Mexican children. Importantly, we found new 
associations between MG4R and fasting glucose levels, and between NPC1 and fasting insulin levels. 
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Background 

Obesity and associated comorbidities (such as cardiovas- 
cular diseases, type 2 diabetes, musculoskeletal disorders 
and some cancers) represent a major public health prob- 
lem that, in recent years, has reached epidemic propor- 
tions. In Mexico, obesity is the most common nutritional 
disorder of children, and the prevalence of this disease 
has increased alarmingly over the past decade: recent 
data from ENSA 2006 ('Encuesta National de Salud', i.e. 
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The National Health Census) reported an increase in 
overweight and obesity from 18% in 1999 to 26% in 2006, 
in young Mexican children (between 5 and 11 years old) 
[1]. Overweight in childhood or adolescence is a risk fac- 
tor for overall mortality in adulthood [2]. 

The prevalence of obesity varies across populations 
and it is noteworthy that the Mexican population has 
been disproportionally affected. Although environmental 
changes (linked to the 'westernization' of ways of life) 
can explain the increase in prevalence of obesity on a 
global level, individual variation of body mass index 
(BMI) persists in the same environment, and the heri- 
tability of this trait is very high. Indeed, heritability esti- 
mates of BMI range between 40% and 70% according to 
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studies [3]. The high genetic susceptibility to obesity 
(and both insulin resistance and type 2 diabetes) in the 
Mexican population may be ascribed to the American 
Indian heritage [4]. 

To date, genome wide association studies (GWAS) and 
linkage studies, mostly performed in European adult 
populations, have identified more than 50 loci associated 
with obesity or BMI [5-10]. However, replication attempts 
have yielded inconsistent outcomes [11-13]. In particular, 
replication in other populations is not obvious, and several 
factors (such as ethnic differences in linkage disequilib- 
rium patterns, ethnic-specific associations, gene x environ- 
ment interactions) may puzzle the picture. Furthermore, 
very few studies have been performed in children present- 
ing with severe obesity in non-European populations. 

In the present study, we aimed to perform a follow-up 
replication study including sixEuropean obesity-associated 
genetic variants, in Mexican children. 

Methods 

Study participants 

In the present study, we analyzed 1,685 children (aged 6 
to 12 years) of Mexican origin, from five different states 
of Mexico (San Luis Potosi, Queretaro, Tijuana, Guana- 
juato and Mexico city), who were randomly selected and 
invited to participate in a cross-sectional study between 
2007 and 2011 from public and private schools. Child 
assent was obtained and parents provided written 
informed consent. We collected data from the children 
and parent or legal guardian per child by direct 
questioning. 

Participants were scheduled for clinical laboratory 
evaluation following a 12 h overnight fasting. Blood sam- 
ples were drawn to assess levels of fasting glucose, fast- 
ing insulin, total cholesterol and triglycerides levels. 
Biochemical variables were measured using an ILab 350 
Clinical Chemistry System (Instrumentation Laboratory 
IL). Weight was measured with a digital scale (Seca) and 
height was measured with a portable stadiometer (Seca 
225). BMI was calculated and classified according to the 
'Centers for Disease Control and Prevention 2000' (CDC 
2000) reference [14]. CDC 2000 growth charts are based 
on 5 U.S nationally representative surveys conducted be- 
tween 1963 and 1994, in which Mexican American chil- 
dren were included [14]. According to those growth 
charts, for ages 2 to 20 years, overweight was defined as 
a BMI-for-age between the 85 th and 95 th percentiles, 
while obesity was defined when BMI-for-age was higher 
than the 95 th percentile [14,15]. Insulin resistance was 
defined as: homeostasis model assessment of insulin re- 
sistance (HOMA-IR = [(Fasting glucose ( mg /dL)) (Fasting 
insulin ((lU /mL))]/405) > 3.4 (that is the 90 th percentile of 
HOMA-IR in a population of healthy Mexican children 
[16,17]). 



The study was approved by the local ethics committees 
('Instituto Mexicano del Seguro Social', i.e. the Mexican 
national health service; reference number: 2011-000-079) 
and is in compliance with the Helsinki Declaration. 

SNPs Selection and genotyping 

To achieve a power of at least 50%, with an odds ratio of 
1.25, we only selected SNPs with minor allele frequen- 
cies > 15% in the Mexican population according to the 
HapMap database. These SNPs were identified by GWAS 
or meta-analyses of GWAS in European populations 
[6,7,9,10]. Furthermore, we selected a SNP that was sig- 
nificantly associated with childhood obesity in a French 
population, according to a linkage association study [8]. 
The expected statistical power per SNP (estimated by 
Quanto software) is reported in Additional file 1: Table SI. 

Genomic DNA was isolated from peripheral blood 
white cells using a QIAamp DNA kit (Qiagen), purity 
and integrity was verified by 260/280 nm measurements 
and by electrophoresis in 0.8% agarose gels, stained with 
ethidium bromide. 

All DNA samples were genotyped for the MC4R 
rsl7782313, NEGRI rs2815752, ENPP1 rs7754561, 
NPC1 rsl805081, GNPDA2 rsl0938397 and FTO 
rsl421085 using TaqMan assays on an ABI7900 system, 
following the manufacturer's protocol (Applied Biosys- 
tems). The genotype success rate was at least 98%, and 
no deviation (P>0.01) from Hardy- Weinberg equilib- 
rium was observed in our population. Thirty duplicate 
quality control samples were included and were geno- 
typed with 100% concordance. 

Statistical analysis 

The effect of SNPs on obesity status was assessed using 
a logistic regression model adjusted for age, gender and 
Mexican state (1: San Luis Potosi, 2: Queretaro, 3: 
Tijuana, 4: Guanajuato and 5: Mexico city), under an 
additive model. We also analyzed the effect of SNPs on 
several metabolic quantitative traits (BMI, fasting serum 
insulin, fasting plasma glucose, total cholesterol and 
triglycerides) using linear regressions under an additive 
model adjusted for age, gender and Mexican state. Data 
for fasting serum insulin and triglycerides were logarith- 
mically transformed before statistical analysis. By apply- 
ing Bonferroni correction, a significant p-value has 
been considered when below 1.4 x 10" (0.05/36) and a 
p-value between 0.05 and 1.4 x 10" 3 has been considered 
as nominally significant. All statistical analyses were per- 
formed using the SPSS software (version 14.0). 

Results 

From an initial sample of 1,685 children (between 6 and 
12 years old), we extracted 949 lean and 514 obese chil- 
dren (see clinical characteristics in Table 1) for the case- 



Mejia-Bem'tez et ah BMC Medical Genetics 2013, 14:21 
http://www.biomedcentral.eom/1 471 -2350/1 4/2 1 



Page 3 of 6 



Table 1 Clinical characteristics of lean and obese Mexican 
children 



Characteristics 


Lean children 
BMI <85pc 
W = 949 


Obese children 
BMI >95pc 
« = 514 


P* 


Female (%) 


48 


.4 


43.6 


0.04 


Age (Years) 


9.5: 


1 1.9 


9.5 ± 1 .8 


0.383 


Anthropometric 










Weight (kg) 


32.6: 


±7.8 


49.7 ±13.4 


<0.001 


Height (cm) 


1 34.3 : 


±11.8 


139.6 ±11.7 


<0.001 


BMI (kg/m 2 ) 


17.7: 


±1.4 


25.0 ±3.4 


<0.001 


BMI (percentile) 


65.8: 


±11.2 


97.6 ± 1 .4 


<0.001 


Metabolic factors 










Fasting Insulin (ull/mL) 


8.6: 


±5.9 


14.1 ±9.4 


<0.001 


Fasting glucose (mmol/L) 


4./: 


±0.6 


4.9 ±0.6 


0.001 1 


Total Cholesterol (mg/dL) 


151.8: 


±29.6 


160.1 ±35.1 


<0.001 


Triglycerides (mg/dL) 


88.5: 


±38.2 


126.6 ±67.6 


<0.001 



Data are means ± standard deviation. 

*P-value according to Student's t test comparing lean and obese children. 
BMI, body mass index; pc, percentile. 



control study, and 1,171 nonobese children (with BMI < 
95 th percentile) for the study of metabolic quantitative 
traits including BMI, fasting serum insulin, fasting 
plasma glucose, total cholesterol and triglycerides. As 
expected, obese children exhibited higher values of insu- 
lin and cholesterol levels (P < 0.001) compared to lean 
children. Insulin resistance was found in 45% of obese 
children. 

In both studies, we genotyped six SNPs that are 
known to be associated with risk of obesity in European 
populations: rsl7782313 near MC4R [6]; rs2815752 near 
NEGRI [7]; rs7754561 near ENPP1 [8]; rsl805081 in 
exon 6 of NPC1 [9]; rsl0938397 near GNPDA2 [7] and 
rsl421085 in intron 1 of FTO [10]. 

We identified a significant contribution of GNPAD2 
rsl0938397 to risk of obesity (G risk allele [in European 

populations]; Odds ratio [95 o /o confidence interval] [OR] = 

l-30[i.ii;i. 5 3j; P = 1.34 x 10~ 3 ; Table 2). Furthermore, we 
found nominal association between risk of obesity 
and both ENPP1 rs7754561 (G risk allele; OR = 



0.84 



; P = 0.020; Table 2) and MC4R rsl7782313 (C 



[0.72;0.97] 

risk allele; OR = 1.40[i.o 6i i.83]; P = 0.016; Table 2). This last 
MC4R SNP also showed a nominal effect on BMI (C risk al- 
lele; P [stan dard error] = 0.41 [016] kg/m 2 ; P = 0.012; Table 3). 
Moreover, we identified a nominal association between 
NEGRI rs2815752 and BMI (A risk allele; (3 = 



Table 2 Association between SNPs MC4R rs17782313, NEGRI rs2815752, ENPP1 rs7754561, A/PC 7 rs1805081, GNPDA2 


rs10938397 and FTO rs1421085, and obesity in Mexican children 






Locus SNP RA* 


Cases/ 


N 


RAF 


OR P 




Controls 




(%) 


[95% CI] 




Lean children 


949 


74.4 




NEGRI rs2815752 A 












Obese children 


514 


77.6 


1.18 [0.99; 1.41] 0.068 




Lean children 


949 


74.6 




NPC1 rs1 805081 A 












Obese children 


514 


74.2 


1 .02 [0.85;1 .21] 0.845 




Lean children 


949 


48.6 




ENPP1 rs7754561 G 












Obese children 


514 


44.2 


0.84 [0.72;0.97] 0.020 




Lean children 


949 


7.3 




MCR4 rs1 778231 3 C 












Obese children 


514 


9.8 


1.40 [1.06;1.83] 0.016 




Lean children 


949 


32.2 




GNPDA2 rs 10938397 G 












Obese children 


514 


38.0 


1.30 [1.11 ;1. 53] 1.34X10" 3 




Lean children 


949 


19.1 




FTO rs 142 1085 C 












Obese children 


514 


21.1 


1.13 [0.93;1 .38] 0.228 



*We noted the risk-alleles reported in European populations. 

RA, risk-allele; RAF, risk-allele frequency; OR, odds ratio; CI, confidence interval; P, P-Value. 
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Table 3 Association between SNPs MC4R rs17782313, NEGRI rs2815752, ENPP1 rs7754561, NPC1 rs1805081, GNPDA2 



rsl 0938397 and FTO rs1421085, and metabolic quantitative traits in 1,171 nonobese Mexican children 


Locus SNP RA* 


0(SE) 


BMI 

P 


Fasting serum insulin 


Fasting plasma glucose 


Total cholesterol 


Triglycerides 


0(SE) 


P 


0(SE) 


P 


P(SE) 


P 


0(SE) 


P 


NEGRI rs2815752 A 


0.24 


0.019 


0.01 


0.611 


0.07 


0.353 


2.14 


0.139 


0.01 


0.425 




(0.10) 




(0.03) 




(0.07) 




(144) 




(0.02) 




NPC1 rsl 805081 A 


-0.16 


0.126 


-0.10 


1.21 x10 -3 


-0.08 


0.278 


2.22 


0.144 


-0.003 


0.884 




(0.11) 




(0.03) 




(0.07) 




(1.52) 




(0.02) 




ENPP1 rs7754561 G 


-0.11 


0.209 


0.004 


0.873 


0.01 


0.905 


2.20 


0.087 


0.004 


0.827 




(0.09) 




(0.03) 




(0.06) 




(1.29) 




(0.02) 




MCR4 rsl 77823 13 C 


0.41 


0.012 


0.04 


0.362 


0.36 


1.47 xlO 3 


2.71 


0.247 


0.02 


0.410 




(0.16) 




(0.05) 




(0.11) 




(2.34) 




(0.03) 




GNPDA2 rsl 0938397 G 


0.10 


0.288 


0.04 


0.126 


0.10 


0.132 


1.90 


0.176 


0.03 


0.076 




(0.10) 




(0.03) 




(0.07) 




(140) 




(0.02) 




FTO rs 142 1085 C 


-0.17 


0.171 


-0.06 


0.064 


-0.14 


0.111 


3.70 


0.032 


0.006 


0.778 




(0.12) 




(0.04) 




(0.09) 




(1.72) 




(0.02) 





Data for fasting serum insulin and triglycerides were logarithmically transformed before statistical analysis. Units are: BMI in kg/m , fasting serum insulin in uU/mL, 
fasting plasma glucose in mmol/L, total cholesterol in mg/dL and triglycerides in mg/dL. 
*We noted the risk-alleles reported in European populations. 



RA, risk-allele; SE, standard error; 0, effect-size; P, P-value; BMI, body mass index. 

0.24(0.10] kg/m 2 ; P = 0.019; Table 3). Of note, we did 
not find any significant association between risk of 
obesity (or BMI variation) and either NPC1 
rs2815752 or FTO rsl421085 (Tables 2 and Table 3). 

In our study of quantitative traits, we found a signifi- 
cant association between NPC1 rsl 805081 and fasting 
serum insulin levels (A risk allele; (3 = -0.10[o.o3] uU/mL; 
P= 1.21 x 10~ 3 ; Table 3). We also identified an effect of 
the obesity risk allele of MC4R rsl7782313 on increased 
fasting plasma glucose levels (C risk allele; (3 = 0.36[o.n] 
mmol/L; P= 1.47 x 10' 3 ; Table 3). A trend of association 
was found between FTO rsl421085 and both total chol- 
esterol levels (C risk allele; (3 = 3.70[i. 72 ] mg/dL; P = 0.032; 
Table 3) and fasting serum insulin levels (C risk allele; 
P = -0.06 [0 .o4] uU/mL; P = 0.064; Table 3). 

Discussion 

From our study based on obese and nonobese Mexican 
children, we found a significant contribution of the 
minor allele of GNPDA2 rsl0938397 to increased risk of 
obesity. The association between GNPDA2 and BMI had 
firstly been identified in a meta-analysis of several Euro- 
pean GWASs performed by the GIANT consortium [7]. 
Subsequently, the association signal has been confirmed 
in other populations: in adults from East-Asia [18,19] 
and in Chinese children [20]. To our knowledge, the 
present study is the first confirmation of the GNPDA2 
association signal with risk of obesity in the Mexican 
population. 

Furthermore, we found nominal associations between 
risk of obesity or BMI and the following SNPs: ENPP1 
rs7754561, MC4R rsl7782313 and NEGRI rs2815752. 



Of note, the risk alleles for obesity or increased BMI 
were the same between Europeans and Mexicans, except 
for the risk allele of ENPP1 rs7754561 (in Europeans) 
that showed a protective effect in Mexicans. Although 
the association between both MC4R and NEGRI and 
risk of obesity has been confirmed in a plethora of stud- 
ies and populations, the ENPP1 association signal with 
obesity is more controversial [21]. To our knowledge, no 
other studies demonstrate a protective role of the risk al- 
lele of ENPP1 rs7754561. Recently, it has been shown 
that ENPP1 overexpression in human adipocyte cell lines 
resulted in defective adipocyte maturation [22]. If con- 
firmed in other Mexican populations, the protective ef- 
fect of the ENPP1 variant may be due to a loss-of- 
function of the protein. 

We did not find any significant contribution of either 
NPC1 rs2815752 or FTO rsl421085 to obesity or BMI 
variation. Therefore, the FTO association signal with 
obesity which was found in Mexican adults by 
Villalobos-Comparan and colleagues [23], was not con- 
firmed in children. Thus, the obesogenic effect of FTO 
would occur later in the Mexican population than in the 
Europeans [9], but as we lack some statistical power, 
additional genetic studies on Mexican children would be 
needed. 

Importantly, we identified a significant effect of the 
risk allele of NPC1 rsl805081 on decreased fasting 
serum insulin levels. NPC1 encodes Niemann Pick dis- 
ease type CI protein that mediates intracellular choles- 
terol trafficking via binding of cholesterol to its N- 
terminal domain. Very recently, Jelinek and colleagues 
have reported that Npcl haploinsufficiency developed 
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abnormal metabolic features (including hyperinsuline- 
mia) and increased susceptibility to weight gain in mice 
[24]. We could suspect that the variant may have a gain- 
of-function effect in Mexican children. Of note, we did 
not report any association between the NPC1 variant 
and obesity or BMI variation, and this association signal 
remains quite controversial according to population 
studies [20,25,26]. 

Furthermore, we identified an effect of the risk allele 
of MC4R variant and increased fasting plasma glucose 
levels. To our knowledge, no previous study has shown 
this association. However, two recent large meta- 
analyses of GWAS identified a significant association 
between the MC4R locus and type 2 diabetes risk, in 
European and Asian populations [27,28]. Of note, our 
present association with fasting plasma glucose levels 
remained significant after adjustment for BMI (data not 
shown). Altogether, these findings would suggest a 
potential effect of MC4R polymorphisms on decreased 
pancreatic beta-cell function. 

Conclusion 

In summary, in our sample of Mexican children, we repli- 
cated four European obesity-related genes (GNPDA2, and 
nominally, NEGRI, MC4R and ENPP1) in the same direc- 
tion of effect as previous findings, except for ENPP1. Inter- 
estingly, we found two novel association signals: between 
NPC1 variant and fasting serum insulin levels, and 
between MC4R variant and fasting plasma glucose levels. 
These findings should deserve some confirmatory studies 
in other Mexican populations. 

It is noteworthy that these well established genetic 
associations with obesity explain very little of the genetic 
risk for pediatric phenotype, suggesting the existence of 
additional loci whose number and effect size remain un- 
known, which guarantee intense additional investigations 
in the near future. 

Additional file 



Additional file 1: Table SI. Expected statistical power for reaching an 
odds ratio of 1.25 in the present obesity case-control study (N casES = 5 1 4; 
Wco„ fro / 5 = 949). 
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